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Abstract. We construct a hadronic kinetic model which describes production of strange particles in ul- 
trarelativistic nuclear collisions in the energy domain of SPS. We test this model on description of the 
sharp peak in the excitation function of multiplicity ratio K + /ir + and demonstrate that hadronic model 
reproduces these data rather well. The model thus must be tested on other types of data in order to verify 
the hypothesis that deconfinement sets in at lowest SPS energies. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 25.75.Dw Particle and resonance production 



1 Introduction 

The central issue in studies with ultrarelativistic heavy 
■ ion collisions is to map the phase diagram of strongly 
interacting matter. From lattice QCD and numerous in- 
vestigations with effective models we know that at high 
enough energy density hadronic matter melts to a phase 
where quarks and gluons are the relevant degrees of free- 
dom. Data on jet quenching at RHIC energies indicate 
\ very clearly that such a deconfined phase has been pro- 
duced there [l]|2]|3]|4J . This does not. however, answer the 
question where is the threshold for deconfinement. 

In order to find the minimum collision energy at which 
the hadronic description of collision dynamics turns inade- 
quate, a natural choice is to study various excitation func- 
tions, i.e. dependences of quantities on the collision energy. 
The expectation here is that a change of the quality of 
the system would demonstrate itself as a non-monotonic 
dependence of some excitation function. A set of interest- 
ing excitation functions has been observed in the energy 
scan at the SPS. These include: i) a sharp peak at pro- 
jectile energy 30 AGeV of the ratio of positive kaon to 
pion multiplicities ("the horn"); ii) a change in the slope 
of number of produced pions per participating nucleon 
("the kink"); and iii) a plateau in the excitation function 
of kaon mean p t which is otherwise increasing function of 
the collision energy ("the step") In addition to this, 

event-by-event fluctuation of the K/tt multiplicity ratio 
grows when the collision energy is lowered and reaches 
maximum at lowest SPS energies, though we have to note 
that there are no data at lower energies from AGS. 

These observed features possibly indicate qualitative 
changes of the system. A possible explanation is that they 
are caused by the onset of deconfinement. Such a positive 
reasoning is insufficient, however, for claims of this kind 



of discovery. In fact, one would have to demonstrate that 
any purely hadronic description of the data fails. In this 
paper we embark on such a programme and address one 
of the excitation functions: the horn. 

Currently, the data have been addressed in many dif- 
ferent approaches. Transport codes generally fail in repro- 
ducing some of the observed multiplicities as functions of 
collision energy and thus fail in ratios 7 . A three-fluid hy- 
drodynamic simulation does not correctly reproduce the 
multiplicity of negative pions [Hj- Statistical hadronisa- 
tion model leads to a broad peak of the K + /n + ratio 
as a function of collision energy which is put in connec- 
tion with transition from baryon-dominated to meson- 
dominated entropy The peak is much broader than 
observed, however. Data are better reproduced in a sta- 
tistical hadronisation fit in which strange species are sup- 
pressed with respect to chemical equilibrium by a stran- 
geness suppression factor |10j . A question appears then, 
what is the mechanism leading to the particular value of 
strangeness suppression factor? 

The data are successfully interpreted in framework of 
the so-called Statistical Model of the Early Stage (SMES) 
The model predicts first order phase transition and 
the existence of mixed phase, which sets in about the 
beam energy of 30 ^4GeV. Another important ingredient 
of the model is the assumption of an immediate chemical 
equilibration of the primordial production of quanta. A 
non-equilibrium kinetic calculation reproducing the data 
which also includes first order phase transition has been 
proposed recently ^2] ■ 

We shall construct hadronic non-equilibrium kinetic 
model and try to reproduce the peak in the ratio of mul- 
tiplicities K + /tt + . It will be shown that these data can 
be interpreted in framework of such a model and there- 
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Fig. 1. (colour online) If energy available for production of strange particles increases with collision energy while the lifetime 
of fireball decreases with increasing collision energy, these two excitation functions might combine into a peak in strangeness- 
to-entropy ratio as a function of collision energy. The ratio of strangeness to entropy is experimentally accessible through the 
ratio of positive kaon to pion multiplicity. 



fore hadronic description must be tested on other types of 
data until it can be safely ruled out. 



2 Strangeness production: cartoon description 

Strange particles are produced in nucleon-nucleon reac- 
tions, though in smaller fraction to non-strange ones than 
in nuclear collisions. We shall assume that the correspond- 
ing surplus of strange particles is produced in secondary 
reactions of hadrons generated in nuclear collisions. Then, 
two aspects are important for strangeness production: avail- 
able energy and time. 

The available energy is microscopically represented 
through energy density. At typical momenta of particles 
in the considered environment, many secondary reactions 
run close to the threshold of kaon and hyperon produc- 
tion. In such a case, production rates depend strongly 
on incoming momenta and densities of reactants, which 
both depend on temperature (if it is defined) . Higher en- 
ergy density (and temperature) implies higher production 
rates. 

Since strangeness is under-represented with respect to 
chemical equilibrium in nucleon-nucleon collisions, its 
amount in the system will increase with time until it sat- 
urates. Thus longer total lifespan of the system would 
also mean larger relative amount of strange particles. This 
picture becomes slightly more complicated in an expand- 
ing and cooling system, but the simple assumption that 
the amount of strange particles grows with time is good 
enough for the cartoon-like explanation offered in this sec- 
tion. 

Now we can imagine the following explanation of the 
K + /tt + horn (Fig.^l: the energy density within the sys- 
tem which is available for production of strange particles 
is an increasing function of the collision energy. We could 
also assume that the total lifespan of the fireball decreases 
as the collision energy grows. These two features might be 
combined in such a way, that the strangeness-to-entropy 
ratio — which is basically measured by the K + /ir + ratio — 
shows a sharp peak as it is indeed observed. In the next 
Section we shall describe this model in more technical 
terms. 




Fig. 2. (colour online) Evolution scenarios and freeze-out hy- 
persurfaces in a space-time diagram for fireballs created at 
lower collision energy with strong stopping and re-expansion 
(solid line) , and higher collision energy in nuclear transparency 
regime (dashed line) . Longer lifetime of a fireball does not nec- 
essarily imply larger extent in longitudinal direction. 



The assumed dependence of the total lifespan on col- 
lision energy may result from an interplay between stop- 
ping and subsequent accelerated expansion. At low col- 
lision energy we expect strong slowing down of the inci- 
dent nucleons. Then the system needs time in order to 
build up expansion from pressure and to expand up to 
the breakup density. At higher collision energy stopping 
is weaker and the incident nucleons continue longitudinal 
movement. Thus expansion is present from the first mo- 
ments of fireball evolution and it takes shorter time to 
get to the critical breakup density. Let us also stress that 
the assumption of shorter lifespan at higher collision en- 
ergy does not clearly contradict any existing data. It does 
not imply larger longitudinal extent of the fireball and a 
larger freeze out volume (Fig. EJ. Hence, it also does not 
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necessarily lead to higher multiplicity at lower collision 
energy. In fact, the correspondence between lifetime and 
the longitudinal size is model-dependent and one usually 
thinks about Bjorken boost-invariant expansion when 
making such a connection, but this is not the generally ap- 
plicable hydrodynamic solution. 



3 The kinetic model 

We shall construct a kinetic model for the production of 
strange particles. Such a calculation requires knowledge 
of the local densities of reacting species. This is naturally 
provided in so-called hydro-kinetic or flavour kinetic mod- 
els |14U15lllfilll7lll8lll9| . It is known, however, that most 
hydrodynamic and kinetic models have problems with re- 
producing the space-time extent of the fireball measured 
by femtoscopy. Therefore, we shall not use such a model in 
order to calculate the space-time evolution of the fireball. 
Instead, we will parametrise the evolution of the density. 
In this way, we do not couple our approach directly to 
the microscopic structure of the matter or to its equation 
of state. On the other hand, we can choose many differ- 
ent parametrisations of the fireball evolution and check 
which of them might lead to results consistent with data. 
Our ansatzes for these parametrisation will be constructed 
with one eye on femtoscopy data and hadronic spectra. 

As we are now only interested in ratios of multiplicities 
it will be sufficient to calculate the densities of individual 
species in the kinetic approach. For densities of kaons with 
positive strangeness we derive 



dpK 
dr 



c^Nk 
dr V 



N K I dV 
~V~ V ~dr 



1 dN K 
V dr 



(1) 



Notice that the first term on the right hand side 1/V dV/ dr 
actually includes the expansion rate which will follow from 
our assumption for density evolution. The second term in- 
cludes the rate of change of the number of kaons. It can 
be divided into two contributions: production rate and an- 
nihilation rate. The former stands for all processes which 
produce a kaon while the latter includes processes which 
destroy kaons. They are determined from cross-sections 
of these processes and evolved densities. Thus the master 
equation for kaon density reads 



dpK 
dr 



Pk -77 



I dV 
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3.1 Production and annihilation 

Densities of K + , K°, K* + , and K*° are evolved accord- 
ing to eq. (J2J. Here we always assume that all particles 
keep their vacuum properties. The gain term has in our 
treatment basically two types of contributions 
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where the sum goes over two-to-N processes leading to 
production of kaons and the second term is for K* decay 
into K. The term in angular brackets is cross-section of 
a process multiplied with relative velocity of the reacting 
species and averaged over the distribution of relative ve- 
locities. For the sake of this averaging we always assume 
thermal distribution of velocities. The annihilation term 
is obtained in a similar way 
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In a real calculation, it is impossible to include all pro- 
cesses that create or destroy a kaon. Only the most im- 
portant ones are included here: 

— Associated production of kaon and a hyperon in reac- 
tions of ttN «-» KY and irA «-» KY. In order to keep 
the detailed balance, these reactions are included in 
both directions, i.e., creating and annihilating kaons. 

— Meson-meson reactions of 7T7t, np, and pp are also in- 
cluded in both directions. 

— K* production from ttK collisions and its decay. 

— Reactions of irY <-> KB. 

— Baryon-baryon reactions which lead to kaons in the 
final state are included only in the gain term. 

For details of parametrisations of all the used cross-sections 
we refer the reader to [23]. Note that no reactions involv- 
ing antibaryons have been taken into account. Since nu- 
clear collisions create a baryon rich environment the error 
due to this simplification is reasonable. It is largest at the 
highest SPS energy (we do not apply this model to higher 
energies) where we estimate it around 10%. 

Other species than kaons with S > are not treated 
explicitly in the kinetic approach. Chemical reactions chang- 
ing their numbers are swift and therefore we can assume 
that non-strange species are kept in chemical equilibrium 
during the whole evolution. For species with S < — which 
include K~ and K° as well as A, £, B, and fl — we as- 
sume that they are in relative chemical equilibrium, i.e., 
that all the strange quarks are distributed in the system 
according to principle of maximum entropy although their 
total number is given by the number of produced strange 
antiquarks contained in kaons. 

3.2 The ansatz for expansion 

We shall assume that the evolution of fireball densities 
consists basically from two periods: initial accelerating pe- 
riod and later scaling expansion with power-law depen- 
dence of density on time. Since baryon number is a con- 
served quantum number, evolution of baryon density ac- 
tually determines the expansion of the volume. We write 
it as 



Pq( t ) 



PQo(l - 

p'qo 
(t-t )° 



br 2 ) 



2 \S 



T < T s 
T > T„ 



(5) 



where pqq, a, b, p'q Q , to, t s , a and S are parameters which 
can be tuned. They can be put in relation to the total life- 
time, initial density, initial density decrease etc. [23j. The 
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subscript Q stands for any conserved quantum number; 
in our treatment we explicitly take care of baryon number 
and the third component of isospin. This parametrisation 
is motivated by simplicity: at the beginning the simplest 
way of encoding the acceleration of expansion is through 
second order polynomial, while the power law in the end 
is commonly used in analysis of correlation radii. 

The energy density is parametrised in a similar way 



e(r) = 



£ 'o 

(t-t >/s 



- br 2 ) : t < t s 

■■ T >T S 



(6) 



so the parameter S represents a simple way of choosing 
the equation of state. 

By exploring a range of parameters we can tune our 
parametrisation between the well known Bjorken 13. and 
Landau |21| hydrodynamic solutions of fireball expansion. 



3.3 Initial and final conditions 

Since kaons are also produced in nucleon-nucleon interac- 
tions, there must be some initial kaon density due to pri- 
mordial kaon production in collisions of incident nucleons. 
It is estimated |20| from a compilation of kaon production 
data in nucleon-nucleon collisions |22| . 

Then, species containing strange quarks must balance 
the strange antiquarks such that the total strangeness van- 
ishes. Among themselves they are in relative equilibrium. 

The most important choice is that of the total life- 
time of the fireball and the initial energy density. This, 
together with the final state densities, basically determines 
the whole evolution scenario. 

The energy density and number densities in the final 
state, i.e., at the end of parametrisations JSJ and J§J, are 
chosen so that they correspond to the measured values. 
The values from experiment were calculated from results 
of chemical freeze out fit with a statistical model 

Thus the construction of our model guarantees that 
we end up with the correct final state values of energy 
density, baryon density, and the density of third compo- 
nent of isospin. It still does not guarantee, however, the 
correct K + abundance, as kaons are produced kinetically. 
But once we get the right amount of kaons, together with 
the densities this implies the correct value of temperature 
and chemical potentials. Consequently, the whole chemical 
composition is correct in such a case. 



3.4 Kaon and A production 

From what has been said it should be now easy to under- 
stand how the production of kaons, antikaons, and lamb- 
das is steered by different knobs. This is important as we 
shall fit the ratios of their multiplicities to pion multiplic- 
ity. The amount of kaons is mainly given by the lifetime 
of the whole system and also by the initial energy den- 
sity. The strange quarks are distributed — in a simplified 
picture for the sake of this explanation — among antikaons 
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Fig. 3. (colour online) The ratios of multiplicities K + /ir + 
(left panel), K~ /tt~ (middle), and A/it (right) calculated in 
scenarios with various initial energy densities and total life- 
times of the fireball for Au+Au collisions at projectile energy 
11.6 AGeV. The plots show dependence of the ratios on the 
total lifetime of the fireball. Different curves correspond to dif- 
ferent initial energy densities. Horizontal lines indicate the la 
intervals around the measured values. 



and hyperons. The distribution of these quarks among an- 
tikaons and hyperons is given in relative chemical equilib- 
rium by the temperature and baryonic chemical potential. 

In summary: lifetime determines the amount of kaons, 
and subsequently temperature specifies how many antika- 
ons and hyperons there will be. 



4 Results 

We have run our hadronic kinetic model for various sets 
of model parameters for Au+Au collisions at projectile 
energy 11.6 ^4GeV (highest AGS energy) and then Pb+Pb 
collisions at 30, 40, 80, and 158 ^IGeV (SPS), so that we 
safely cover the region where the K/ir horn appears. For 
each value of initial energy density and the total lifetime 
we evolve kaon densities. At the end, we add the feed-down 
to kaons and pions due to resonance decays. (Details of 
this procedure can be found in |2()j.) 

Some selected results are shown in Figs. |3| 0J and [5] 
In general, we observe that the resulting ratios mainly 
depend on the total lifetime and less so on the initial en- 
ergy density. It seems that the data values can be repro- 
duced with rather reasonable lifetimes, though somewhat 
longer than one would infer from analysis of correlation 
radii based on Bjorken hydrodynamic solution. 

In order to make our comparison with data more quan- 
titative, we calculate for each set of results 



x 2 = 



E 



(7) 



where ti, di, and Oi are the calculated and measured val- 
ues and experimental error, respectively. For all calculated 
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Fig. 4. (colour online) Same as Fig.[3]but for Pb+Pb collisions 
at projectile energy 30 ^4GeV. 
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Fig. 5. (colour online) Same as Fig.[3]but for Pb+Pb collisions 
at projectile energy 158 AGeV. 



scenarios the values of \ 2 are summarised in Fig-EI We ob- 
serve that the maximum of K + /ir + ratio of multiplicities 
as a function of collision energy is basically translated into 
a maximum lifetime of the fireball for the best reproduc- 
tion of data, although such a conclusion is not statistically 
significant. One could speculate if it is connected with long 
total lifetime due to a soft point in the equation of state 
which would appear close to the phase transition. 

On the other hand, based on comparison with data one 
neither can exclude that the total lifetime of the fireball 
does not increase and rather decreases when the collision 
energy grows. In order to demonstrate this we choose for 
each collision energy one set of parameters, and show re- 
sults calculated with those parameters compared to data 
in Fig. The lifetimes and initial energy densities from 
the used parameter sets are collected in Tabled 

From the table we also see that the final state tem- 
perature which we obtain does not differ much from the 
result of chemical freeze-out analysis. Hence, we conclude 
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Fig. 6. (colour online) Quantity \ 2 defined in eq. Q evaluated 
for results from all simulations. Plots show the dependence on 
total lifetime and different curves correspond to different initial 
energy densities. 
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that not only kaon production, but the whole calculated 
chemical composition is correct. 
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Table 1. Initial energy densities and total lifetimes from pa- 
rameter sets which were used in calculations leading to results 
shown in Fig. In the lower portion of the table, T/ is the 
final state temperature obtained in our simulations and T the 
temperature from the analysis of chemical freeze-out |10) . 
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5 Conclusions 

The actual aim in search for the onset of deconfinement 
is to falsify hadronic description of data. This has not 
been accomplished here. Our hadronic kinetic model was 
able to describe the excitation function of the multiplicity 
ratios K + /n + , K~ /ir~ , and A/%. It requires, certainly at 
least for beam energies above 30 AGeV, that the lifetime 
of fireball decreases as the energy of collisions grows. 

So far, we only checked our model against one type of 
data. In order to verify or falsify it, it will have to be tested 
on other interesting excitation functions mentioned in the 
Introduction as well as the "standard" data like, trans- 
verse momentum spectra and correlation radii. Since we 
make statements about the whole evolution of the fireball, 
it seems important to calculate our prediction for dilepton 
spectra which are produced during the whole evolution. 

Finally, so far we have put aside the question of how 
the evolution which we parametrised can result from mi- 
croscopic structure of the matter. If the model passes all 
data tests, this question will have to be addressed. 
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